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Abstract 

Adding  distributed  generation  (DG)  is  a  desirable  strategy  for  providing  highly  efficient  and  environmentally  benign  services  for  electric  power, 
heating,  and  cooling.  The  interface  between  a  solid  oxide  fuel  cell  (SOFC),  typical  loads,  and  the  electrical  grid  is  simulated  in  Matlab/Simulink 
and  analyzed  to  assess  the  interactions  between  DG  and  the  electrical  grid.  A  commercial  building  load  profile  is  measured  during  both  steady-state 
and  transient  conditions.  The  load  data  are  combined  with  the  following  models  that  are  designed  to  account  for  physical  features:  a  One-Cycle 
Control  grid-connected  inverter,  a  One-Cycle  Control  active  power  filter,  an  SOFC,  and  capacitor  storage.  High  penetration  of  DG  without  any 
power  filter  increases  the  percentage  of  undesirable  harmonics  provided  by  the  grid,  but  combined  use  of  an  inverter  and  active  power  filter  allows 
the  DG  system  interconnection  to  improve  the  grid  tie-line  flow  by  lowering  total  harmonic  distortion  and  increasing  the  power  factor  to  unity. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  U.S.  Department  of  Energy  (DOE)  projects  that  total  US 
electricity  demand  will  increase  to  5220  billion  kWh  in  2025 
from  3481  billion  kWh  in  2003  [1].  There  are  numerous  eco¬ 
nomic  and  societal  incentives  for  increasing  the  grid  capacity  to 
meet  this  growth  without  major  changes  to  the  existing  transmis¬ 
sion  infrastructure  [2,3].  Adding  distributed  generation  (DG)  is 
a  viable  way  to  increase  grid  capacity  while  providing  additional 
benefits  over  traditional  electricity  generation  and  transmission 
[2]. 

Distributed  generation,  which  is  the  process  of  generating 
electricity  at  or  near  the  point  of  its  use,  can  both  reduce  the 
demand  for  transmission/generation  and  provide  backup  power 
[3].  Benefits  of  keeping  generation  and  load  physically  close 
include  high  power  quality,  reduced  transmission  losses,  reduced 
dependence  on  the  transmission  system,  and  improved  efficiency 
through  combined  heat  and  power  systems.  Furthermore,  DG 
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provided  by  fuel  cells  could  have  lower  environmental  impact 
and  higher  electrical  efficiency  compared  to  traditional  large 
power  plants  [4,5]. 

One  of  the  first  economically  viable  applications  for  fuel 
cell  DG  is  supplying  the  electrical  load  for  commercial  and 
industrial  applications  that  value  reliable  and  high  quality  power 
[4,5].  Before  distributed  generation  can  be  widely  embraced  as 
a  future  power  technology,  it  must  be  demonstrated  as  beneficial 
to  both  the  user  and  the  connected  power  system.  Several  pre¬ 
vious  research  efforts  have  shown  control  and  load-following 
ability  of  DG  systems  [3,6,7].  This  paper  evaluates  a  sam¬ 
ple  DG-grid  interconnection  on  the  basis  of  experimentally 
measured  steady- state  and  dynamic  commercial  building  load 
data.  The  performance  of  both  the  DG  and  the  grid  is  evalu¬ 
ated  by  investigating  the  effects  of  real  building  load  demand 
on  the  grid  tie-line  current  at  the  point  of  common  coupling 
(PCC). 

The  individual  components  of  the  model  are  first  described 
in  Section  2.  The  methodology  for  load  data  experiment  is 
described  along  with  the  relevant  results.  The  DG  components 
are  a  set  of  dynamic  models  created  within  Matlab/Simulink® 
that  are  designed  to  account  for  the  physical  features  of  a  solid 
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oxide  fuel  cell  (SOFC),  One-Cycle  Control  (OCC)-based  grid- 
connected  inverter  (GCI),  and  active  power  filter  (APF).  In 
Section  3,  the  response  of  each  of  these  models  is  shown  and 
compared  to  previous  experimental  results.  Because  directly 
investigating  the  relationship  between  the  SOFC  and  the  grid 
obfuscates  the  role  of  the  inverter,  Section  4  describes  the 
inverter-grid  relationship  directly  under  steady- state  conditions 
both  with  and  without  an  active  power  filter.  After  steady- state 
operation  is  explored,  the  dynamic  SOFC  model  and  dynamic 
load  data  are  added  to  the  model  in  Section  5  to  investigate  grid 
impacts  under  different  conditions. 

2.  Development  of  load  data  and  models  of  SOFC, 
OCC-GCI,  and  OCC-APF 

2.7.  Load  data 
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Data  from  a  typical  commercial  Southern  California  office 
building  is  acquired  to  provide  insight  into  the  steady- state  and 
transient  power  demand  conditions  that  distributed  generation 
might  be  required  to  meet.  The  3 -phase  power  entering  the 
building  is  monitored  with  a  Nexus  model  1270  power  qual¬ 
ity  monitor,  which  can  measure  both  the  voltage  and  current 
waveforms  with  user-defined  sampling  rates. 

In  order  to  capture  detailed  load  transient  data  two  sampling 
rates  are  employed.  Low-frequency  sampling  was  used  to  cap¬ 
ture  long-term  data,  while  high  frequency  sampling  was  used 
only  to  capture  the  worst-case  load  transient.  Both  the  daily 
power  data  and  the  worst-case  transient  dynamic  data  are  pre¬ 
sented  in  Fig.  1.  There  are  two  levels  of  power  consumption: 
in-use  during  the  workday  and  non-use  at  night  and  weekends. 
As  seen  in  the  close-up  of  Fig.  1,  the  power  demand  does  not 
instantly  jump  from  night  load  to  day  load,  but  is  composed 
of  a  series  of  load  increases  as  each  of  the  several  air  condi¬ 
tioning  motors  turn  on.  The  measured  transient  sample  is  the 
first  and  largest  of  these  start-up  transient  building  demand 
curves,  when  the  building  air  conditioning  system  first  turns 
on  for  the  day.  The  3 -phase  current  and  voltage  waveforms  dur¬ 
ing  this  start-up  are  sampled  at  128  samples/cycle,  or  7680  Hz, 
and  the  resulting  instantaneous  power  calculation  is  shown  in 
Fig.  2. 


Time  of  Day 


Fig.  1.  Diurnal  power  demand  for  5/24/06  with  5  s  sampling  interval  including 
close-up  of  early  AM  transient  with  2  s  sampling  interval. 


Fig.  2.  Dynamic  load  data  measured  during  early  AM  transient  and  sampled  at 
7680  Hz. 

The  instantaneous  power  exhibits  a  high  frequency  oscillation 
because  the  three  fundamental  phases  of  current  average  out  to 
a  constant  value,  but  the  harmonics  do  not. 

2.2.  Solid  oxide  fuel  cell 

The  fuel  cell  type  assessed  in  this  work  is  a  solid  oxide  fuel 
cell.  SOFC  technology  is  being  rapidly  advanced  through  sig¬ 
nificant  public  and  private  investment  [8-10].  The  SOFC  model 
presented  here  is  developed  in  Matlab/Simulink®  according  to 
the  same  methodology  from  [11-14],  but  simplified  for  integra¬ 
tion  with  power  electronics.  This  same  methodology  was  used  to 
develop  models  that  compared  favorably  to  experimental  testing 
of  dynamic  single  cell  transients  [15],  integrated  simple  cycle 
SOFC  systems  [11],  SOFC-MTG  hybrid  systems  [13],  and  PEM 
stationary  fuel  cell  systems  [16]. 

The  SOFC  response  is  fundamentally  limited  by  the  perfor¬ 
mance  of  the  fuel  preprocessor  and  the  amount  of  hydrogen 
present  in  the  anode  compartment  [7,1 1,17].  The  SOFC  system 
model  presented  herein  is  therefore  simplified  by  only  captur¬ 
ing  fuel  transient  effects  within  the  fuel  processor  and  fuel  cell, 
which  limits  SOFC  system  transient  performance.  The  key  equa¬ 
tions  that  capture  fuel  cell  electrochemistry,  hydrogen  transients, 
and  controls  that  are  used  in  the  present  model  are  stated  below. 
The  fuel  cell  voltage  shown  in  (1)— (5)  is  a  common  represen¬ 
tation  found  in  [18].  The  cell  voltage  is  the  Nernst  potential 
(^Nernst)  minus  the  activation  (7s  Act  X  ohmic  (TsohmicX  and  con¬ 
centration  (Ticonc)  polarizations. 
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Hydrogen  transients  within  the  fuel  cell  are  captured  with  the 
dynamic  species  conservation  equation  for  hydrogen  as  in  (6). 


N 


dXH2 

d  t 


—  A^in  AW*H2  A^Consumed 


(6) 


The  amount  of  reacted  hydrogen  is  proportional  to  current  by 
Faraday’s  law  as  in  (7). 


A^ consumed  —  ^ 

nF 


(7) 


During  transients  it  is  essential  that  sufficient  hydrogen  be 
maintained  in  the  fuel  cell  to  sustain  the  fuel  cell  voltage  and 
avoid  damage  caused  by  hydrogen  starvation.  Hydrogen  is  sus¬ 
tained  in  the  anode  compartment  by  controlling  the  fuel  flow  rate 
in  proportion  to  consumed  hydrogen  and  a  desired  utilization. 
Utilization  is  the  amount  of  hydrogen  consumed  compared  to 
the  total  hydrogen  flow  rate  as  follows: 


U  = 


N consumed 

N\  n 


(8) 


From  (7)  the  hydrogen  consumption  is  known  from  the  cur¬ 
rent.  Therefore  the  required  flow  rate  to  maintain  utilization  can 
be  calculated  directly  as  in  (9). 


This  control  strategy  is  called  current-based  fuel  control  and 
is  discussed  in  detail  in  [  1 1  ] .  To  avoid  hydrogen  starvation  during 
flow  delay,  as  explained  below,  i  is  evaluated  from  a  look-up  table 
that  approximates  the  current  demand  from  the  power  demand. 

While  power  demand  can  change  instantaneously,  there  is  a 
finite  time  delay  for  fuel  flow  to  increase.  During  severe  tran¬ 
sients,  this  delay  could  cause  hydrogen  starvation  if  the  current 
drawn  exceeds  the  available  hydrogen.  Limiting  the  fuel  cell 
current  can  thus  prevent  hydrogen  starvation.  The  maximum 
allowable  current,  /max,  is  defined  as  the  current  that  generates  a 
set  maximum  utilization,  Umax,  for  a  given  inlet  hydrogen  flow 


rate,  as  shown  in  (10). 

*max  =  A^in ,  delayed  U max  n  F  (10) 

The  flow  of  hydrogen  to  the  fuel  cell  cannot  increase  instan¬ 
taneously  because  of  delays  in  fuel  delivery,  preparation,  and 
processing.  The  overall  delay  can  include  flow  delays  in  valve 
actuation  and  fuel  desulfurization,  delays  caused  by  the  reformer 
chemical  kinetics  including  those  from  insufficient  heat  trans¬ 
fer  to  the  reformer  and  insufficient  reformer  residence  time,  as 
well  as  mass  flow  delays  due  to  flow  restrictions  in  the  fuel  pro¬ 
cessing  system.  A  reformer  reactor  with  sufficient  heat  transfer 
and  residence  time  is  assumed.  Flow  delays  associated  with  the 
actuators,  the  desulfurizer,  and  the  entire  fuel  processing  system 
are  lumped  into  a  single  fuel  processor  flow  delay.  The  value  of 
the  fuel  flow  actuation  delay  for  a  fuel  cell  system  depends  upon 
the  fuel  processor  design  and  control.  Results  from  [7,19-22] 
indicate  fuel  processor  delay  times  are  generally  on  the  order  of 
seconds.  Thus,  the  current  analyses  first  use  a  hydrogen  delivery 
delay  of  2  s  and  include  a  sensitivity  analysis  to  demonstrate  the 
effects  of  fuel  delay  on  the  response  and  transient  capability  of 
SOFC  systems. 

2.3.  OCC  inverter 

The  inverter  model,  which  is  originally  developed  in  [23]  and 
modeled  in  [24],  has  a  standard  3 -phase  half-bridge  power  stage 
with  an  OCC-based  controller  [25].  Six  switches,  S^p,  •  •  •>  Scn 
control  the  3  phases,  and  inductors  La,  Lb,  and  Lc  connect  the 
half-bridges  to  the  grid.  Vdc  is  the  voltage  at  the  input  DC  bus, 
and  currents  and  ic  are  injected  onto  the  grid.  The  circuit 
is  constructed  in  Matlab/Simulink  in  [24]  with  the  switching 
flow-graph  method  originally  devised  in  [26].  The  circuit  and 
large-signal  model  of  the  power  stage  are  shown  in  Fig.  3. 

The  grid  voltage  and  current  is  sensed  and  processed  by  the 
OCC  controller.  The  complete  360°  3 -phase  voltage  window 
can  be  divided  into  60°  regions,  during  each  of  which  the  duty 
ratio  of  only  two  switches  is  actively  controlled  [23].  The  OCC 
controller  uses  grid  voltage  input  to  select  the  two  active  switches 
and  calculates  the  associated  duty  ratios  Jp  and  dn  based  on  the 


Vm 


Fig.  3.  Power  stage  circuit  diagram  of  OCC  inverter  (a)  and  large-signal  switching  flow-graph  model  of  OCC  inverter  (b)  [24]. 
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Fig.  4.  Power  stage  circuit  diagram  of  OCC-APF  (a)  and  large-signal  switching  flow-graph  model  of  OCC  inverter  (b)  [31]. 


key  control  equation  in  (11)  [23]. 


tion  (14)  in  [28]. 
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Rs  is  the  equivalent  sensing  resistance,  K  is  a  near-constant 
that  limits  output  current,  Vm  is  an  introduced  variable  to  control 
output  power,  zp  and  in  are  selected  injection  currents,  and  Vp 
and  Vn  are  selected  linear  combinations  of  grid  voltages  Va,  Vb, 
and  Vc.  The  key  control  equations  for  K  and  Vm  are  originally 
derived  in  [24]  and  further  developed  in  [27]. 


K  =  1.1  Rs 
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Pref  is  the  desired  output  power  and  Va,rms  is  the  rms  of  input 
grid  voltage  Va. 


2.4.  Active  power  filter 

The  shunt  active  power  filter  (APF)  compensates  current  har¬ 
monics  from  a  non-linear  load,  balances  uneven  phases,  and 
improves  the  power  factor  to  one  [28-30] .  The  OCC-APF  model 
is  originally  derived  in  [31]  by  the  same  control  and  model¬ 
ing  technique  as  the  inverter.  The  APF  power  stage  also  has  6 
switches  Sap, . . .,  Sen  that  are  configured  in  a  3 -phase  half-bridge 
circuit.  V0  is  the  voltage  across  the  capacitor  C,  and  La,  Lb,  and 
Lc  connect  the  three  branches  to  the  grid.  Va,  Vb,  and  Vc  are 
the  grid  voltages  at  the  point  of  contact,  and  /a,  i\>,  and  ic  are  the 
incoming  grid  currents.  The  power  stage  circuit  and  large- signal 
model  are  shown  in  Fig.  4. 

Each  of  the  six  60°  regions  only  has  two  active  switch  pairs, 
and  the  signals  to  each  pair  are  complementary.  The  OCC  con¬ 
troller  uses  the  input  voltages  to  select  which  switches  are  active 
and  generates  the  duty  ratios  dp  and  dn  by  the  key  control  equa- 
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Rs  is  the  sensing  resistance,  Vm  is  a  compensated  feedback  of 
the  voltage  control  loop,  and  ip  and  in  are  selected  input  currents. 


3.  Verification  of  component  models 

Before  they  are  integrated  into  a  DG  system  model,  the  behav¬ 
ior  of  each  of  the  previously  described  physical  models  is  shown 
to  illustrate  their  individual  behavior.  The  inverter  and  APF  mod¬ 
els  have  both  been  previously  built  and  tested  experimentally. 
The  component  model  results  shown  here  match  the  experimen¬ 
tal  results.  Fig.  5  shows  the  response  of  the  inverter  to  a  step 
change  in  power.  The  dynamic  rapidly  adjusts  to  the  power 
increase  at  5  s,  and  the  current  is  sinusoidal  and  in  step  with 
the  voltage  throughout.  Thus  the  inverter  produces  3 -phase  cur¬ 
rent  at  unity  power  factor  and  low  harmonic  distortion,  which  is 
shown  by  experiments  in  [23,24]. 


Fig.  5.  Response  of  inverter  to  step  power  change:  (top  to  bottom)  power  to 
inverter,  phase  A  voltage,  phase  A  current. 
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Fig.  6.  Response  of  shunt  APF  to  a  non-linear  load  for  phase  A:  (top  to  bottom) 
line  voltage,  compensated  line  current,  APF  compensation  currents,  and  load 
current. 


Simulation  results  for  the  APF  model  are  shown  in  Fig.  6. 
A  non-linear  load  is  simulated  to  match  the  experimental  load 
from  [32].  The  load  current  for  phase  A  is  shown  at  the  bottom 
of  Fig.  6,  and  the  line  voltage  for  phase  A  is  shown  at  the  top. 
The  second  graph  shows  the  phase  A  grid  line  current  after  com¬ 
pensation,  and  the  third  graph  shows  the  phase  A  compensation 
current  from  the  shunt  APF.  All  simulation  results  match  with 
the  experiment  from  [32]. 

The  SOFC  model  response  is  presented  in  Fig.  7.  When  the 
power  demand  increases  at  5  s,  the  SOFC  stack  does  not  have 
enough  fuel  to  meet  the  new  demand.  The  utilization  increases 
to  its  maximum  value  of  90%,  which  allows  for  a  slight  increase 
in  output  power,  but  the  SOFC  cannot  match  the  power  demand 
until  7  s  when  the  higher  fuel  rate  can  be  delivered  to  the  anode 
compartment.  This  is  a  direct  result  of  the  2  s  fuel  flow  delay  in 
the  model. 


== 

3-0 

Inverter 

ljnv,a 

rvvv^ 

lline,a  V  Hne,a 

Ijnv.b 

*  \/ 

lline,b  *  lme,b 

ljnv,c  ^ 

> - 

Line,c  V  line,c 

— f  © — 1 

nP 

r  demand 

%-r 

Pgi-id 


Load, a 


load 


}load,b 


jload,c 


Nonlinear 

Load 


Fig.  8.  Block  diagram  of  control  scheme  for  calculation  of  inverter  power 
demand  and  grid  connection. 

4.  Steady-state  infinite  load-following 

The  baseline  case  for  the  grid-connected  inverter  observes  the 
grid-inverter  behavior  independent  of  the  fuel  cell.  The  power 
demand  is  measured  from  the  load,  reduced  by  lOkW  to  avoid 
exportation  back  to  the  grid,  and  input  to  the  inverter.  To  provide 
power  in  the  absence  of  a  prime  mover,  an  infinitely  large  DC 
capacitor  is  assumed  across  the  DC  terminals  of  the  inverter. 
The  grid  current  is  calculated  instantaneously  as  the  difference 
between  building  load  current  and  injected  inverter  current.  The 
block  diagram  is  shown  in  Fig.  8  and  the  steady-state  simulation 
results  are  shown  in  Fig.  9. 

Tie-line  current  is  the  waveform  supplied  by  the  grid  and 
it  contains  large  amounts  of  highly  undesirable  distortion  as 
shown  in  the  bottom  of  Fig.  9.  This  distortion  occurs  because 
the  inverter  follows  the  voltage  reference  and  is  not  designed  to 
supply  the  higher-order  harmonics  that  result  from  a  non-linear 
load.  So  the  inverter  produces  the  majority  of  the  fundamental 
frequency  and  the  grid  supplies  the  harmonics.  Adding  a  shunt 
active  power  filter  as  presented  in  Fig.  10  eliminates  this  problem 
by  providing  compensation  currents  locally  instead  of  from  the 


Fig.  7.  Response  of  SOFC  model  to  step  load  demand  change:  (top)  power  Fig.  9.  Inverter-only  results  for  phase  A:  (top  to  bottom)  line  voltage,  load 
demand  and  SOFC  power  and  (bottom)  fuel  utilization.  current,  inverter  current,  and  line  current. 
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Fig.  10.  Block  diagram  of  control  scheme  with  inverter  and  APF. 


Fig.  12.  Block  diagram  of  DG  system  with  no  load-following  strategy. 


utility.  Additional  benefits  of  adding  an  APF  is  that  it  will  correct 
the  power  factor  of  imported  power  to  1 ,  and  balance  the  phases. 

The  APF  and  inverter  can  be  integrated  into  a  single  unit  with 
One-Cycle  Control,  so  adding  this  power  quality  stage  would  not 
require  a  substantial  increase  in  cost  or  size.  The  line  current  after 
APF  implementation  is  shown  in  Fig.  1 1  and  shows  exemplary 
low  distortion  current  in  the  grid  current. 

5.  Dynamic  grid  behavior  with  a  solid  oxide  fuel  cell 

5.7.  SOFC  for  load-following 

Section  4  showed  that  a  grid-beneficial  scenario  could  be 
implemented  through  the  use  of  an  inverter  and  active  power 
filter  in  the  steady- state.  The  next  concern  is  the  original  goal 
of  observing  grid  behavior  during  a  transient  load  change  with 
a  finite  prime  mover.  The  distributed  energy  resource  used  here 
is  a  SOFC  that  cannot  immediately  follow  the  power  increase 
because  of  the  delay  in  delivering  more  fuel  to  the  fuel  cell  stack. 
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The  block  diagram  of  the  SOFC  incorporation  into  the  system 
is  presented  in  Fig.  12. 

The  inverter  calculates  the  power  reference  and  low-pass  fil¬ 
ters  it  to  protect  the  SOFC  from  the  high-frequency  oscillations 
observed  in  Fig.  2.  The  SOFC  receives  the  power  demand,  and 
the  inverter  injects  the  available  3 -phase  current  onto  the  grid. 
The  deficiency  between  power  delivered  to  the  inverter  and 
power  demanded  is  shown  at  the  top  of  Fig.  13.  The  power  short¬ 
fall  between  load  and  SOFC  is  compensated  by  a  large  increase 
in  grid  tie-line  current  that  creates  a  surge  away  from  the  desired 
grid  power  import.  The  tie-line  current  returns  to  standard  lev¬ 
els  as  the  SOFC  power  rises  to  match  the  load  demand,  but 
the  existence  of  this  spike  indicates  poor  behavior  for  dynamic 
load-following.  To  alleviate  this  load-following  problem,  a  DC 
capacitor  is  added  to  the  inverter. 

5.2.  SOFC  with  capacitors  for  load-following 

Adding  a  large  capacitor  to  the  DC  input  of  the  inverter  aug¬ 
ments  the  load-following  capability.  The  capacitor  is  charged 
during  steady- state  operation  and  discharged  during  the  transient 


Fig.  1 1 .  Inverter  with  APF  results  for  phase  A:  (top  to  bottom)  line  voltage,  load  Fig.  13.  SOFC  with  inverter  and  APF:  (top)  power  demand  and  SOFC  power 

current,  inverter  current,  and  line  current.  and  (bottom)  3-phase  rms  tie-line  current. 
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power  increase  to  supplement  the  SOFC  power.  The  proposed 
control  from  [27]  is  used  to  control  the  voltage  across  a  capac¬ 
itor.  As  the  capacitor  voltage  decreases,  the  difference  between 
reference  voltage  and  actual  voltage,  A  Vc,  is  fed  to  a  PID  con¬ 
troller  that  increases  the  power  demand  to  the  fuel  cell.  When 
the  fuel  cell  is  capable  of  meeting  the  power  demand,  the  sur¬ 
plus  power  will  recharge  the  capacitor.  The  block  diagram  of 
this  system  is  shown  in  Fig.  14. 

The  capacitance  needed  to  completely  supply  the  power 
shortfall  during  the  transient  was  determined  by  simulating  a 
wide  range  of  values  until  the  smallest  value  that  allowed  the 
inverter  power  to  match  the  power  demand  was  found.  This 
value,  which  was  0.96  F,  is  designated  the  full  capacity  capaci¬ 
tor,  and  half  of  this  value,  or  0.48  F  is  the  half  capacity  situation. 
Fig.  15  shows  the  full  capacity  scenario:  the  inverter  power  can 
continuously  match  the  power  demand  (top),  and  the  grid  cur¬ 
rent  is  constant  throughout  the  entire  dynamic  ramping  region 
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Fig.  15.  SOFC  with  inverter,  APF,  and  100%  DC  capacitor:  (top)  power  demand 
and  SOFC  power  and  (bottom)  3-phase  rms  tie-line  current. 


Table  1 

Per-phase  rms  load  current  before  addition  of  DG 


Time  (s) 

A, rms 

7b, rms 

7c, rms 

1.00 

73.08 

57.28 

51.53 

5.00 

92.18 

81.11 

74.69 

8.00 

96.43 

81.41 

76.66 

(bottom).  The  overshoot  in  SOFC  power  is  the  recharging  of  the 
depleted  DC  capacitor. 

An  additional  result  shows  the  value  of  phase  balancing  on 
the  system.  Table  1  shows  the  rms  load  current  for  each  phase 
at  selected  times  before,  during,  and  after  the  transient.  The 
current  magnitude  not  only  increases  during  the  transient,  but  the 
3  phases  are  clearly  unbalanced.  Table  2  shows  the  rms  currents 
at  the  same  times  for  the  system  that  incorporates  APF/Inverter 
DG  along  with  a  DC  capacitor.  The  current  magnitude  remains 
constant  and  the  phases  are  balanced. 

Similar  results  for  a  smaller  half-capacity  DC  capacitor  are  in 
Fig.  1 6.  At  first  the  inverter  power  can  match  the  demand,  but  the 
capacitor  voltage  cannot  drop  too  far  because  it  is  directly  con¬ 
nected  to  the  inverter  terminals.  After  the  capacitor  is  exhausted 
and  before  the  SOFC  can  increase  output,  the  inverter  power 
drops  and  results  in  an  increase  in  tie-line  grid  current.  This 
current  spike  is  less  severe  than  for  the  case  that  is  completely 
without  capacitor  storage,  but  it  does  not  entirely  meet  the  goal 
of  a  constant  tie-line  current. 


Table  2 

Per-phase  rms  line  current  with  APF/inverter  system 


Time  (s) 

7a, rms 

7b, rms 

7c,  rms 

1.00 

8.50 

8.55 

8.53 

5.00 

8.90 

8.95 

8.94 

8.00 

8.36 

8.40 

8.38 
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Fig.  16.  SOFC  with  inverter,  APF,  and  50%  DC  capacitor:  (top)  power  demand 
and  SOFC  power  and  (bottom)  3-phase  rms  tie-line  current. 

Table  3 


Per-phase  rms  line  current  with  APF/inverter  system 


Fuel  delay  (s) 

DC  capacitance  (F) 

1 

0.21 

2 

0.96 

3 

1.51 

4 

2.07 

5 

2.64 

The  fuel  flow  delay  is  the  main  factor  that  determines  the 
rate  of  response  of  the  SOFC  to  a  load  change.  As  this  resulting 
power  shortfall  determines  the  capacitor  size,  the  estimation  of 
fuel  flow  delay  will  correlate  with  the  necessary  capacitance  as 
shown  by  a  sensitivity  analysis  of  fuel  flow  delay.  The  mini¬ 
mum  capacitance  needed  for  fuel  delays  from  1  to  5  s  is  shown 
in  Table  3.  The  SOFC  response  corresponding  to  delays  of  1 
to  3  s  is  also  shown  along  with  the  power  demand  in  Fig.  17. 
The  overshoot  in  SOFC  power  increases  with  fuel  flow  delay 
because  the  larger  associated  capacitor  requires  more  power  for 


Fig.  17.  Response  of  the  SOFC  power  to  load  power  demand  for  fuel  flow  delay 
times  of  1,  2,  and  3  s. 


recharging.  In  all  of  these  cases,  the  charging  and  discharging 
of  the  DC  capacitor  allow  the  inverter  power  to  match  the  load 
demand. 

6.  Conclusions 

This  paper  has  shown  the  development  of  a  model  designed 
to  assess  the  effects  of  a  commercial  building  dynamic  load 
that  incorporates  DG  on  the  electric  power  system  at  the  PCC, 
with  the  goal  of  grid-friendly  behavior.  The  results  show  that 
appropriate  design  of  the  DG  interconnection  system  can  pro¬ 
vide  additional  benefits  for  the  fuel  cell  DG  and  electric  power 
system.  Several  physical  inabilities  of  the  fuel  cell  can  be  allevi¬ 
ated  or  eliminated  with  modest  additions  to  the  inverter  system. 
Some  of  these  inabilities  addressed  here  include:  inability  for 
perfect  load-following,  inability  to  meet  harmonics-based  high 
frequency  power  demand  oscillations,  and  inability  to  provide 
harmonics  compensation.  These  inabilities  are  addressed  at  the 
inverter  through  the  addition  of  an  active  power  filter  and  electric 
energy  storage. 

The  proposed  electric  energy  storage  strategy  includes  a 
large  DC  capacitor,  which  may  be  too  large  for  actual  field 
implementation.  Thus  alternative  storage  methods,  such  as  ultra¬ 
capacitors,  batteries,  and  flywheels,  which  all  have  higher  energy 
density,  should  be  further  explored.  Additional  concerns  include 
a  load-following  ability  and  fuel  flow  delay  attributed  to  the 
SOFC  that  may  be  impractical  for  first  generation  systems.  One 
of  the  goals  of  investigating  DG  systems  is  to  identify  what 
features  next  generation  SOFC  systems  will  need  to  become 
successfully  implemented. 

The  OCC-based  power  electronics  allow  a  system  consisting 
of  a  non-linear  load  and  distributed  generation  to  be  seen  as  a 
simple,  constant  emulated  resistor,  even  during  rapid  dynamic 
transients  of  prime  mover  and/or  load. 

Appendix  A.  Appendix 

The  nomenclature  for  the  SOFC  model  is  as  follows: 

A  cell  surface  area  (m2) 

E  polarization  (V) 

Eo  EMF  at  standard  pressure:  0.92  V  for  hydrogen  and 
oxygen  at  1000  °C 

F  Faraday’s  constant:  96,487  C  mol-1 
i  electrical  current  (A) 

z'o  exchange  current  density  (Am-2) 

il  limiting  current  density  (A  m-2) 

n  number  of  participating  electrons  in  the  reaction 

N  molar  capacity  (kmol) 

N  molar  flow  rate  (kmol  s  - 1 ) 

P  pressure  (kPa) 

R  ideal  gas  constant:  8.314  J  K-1  mol-1 
Rq sr  equivalent  series  resistance  (Q) 

T  temperature  (K) 

U  fuel  utilization 

V  voltage  (V) 

X  species  mole  fraction 
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